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Abstract We review our recent studies demonstrating that the radiation pressure 
problem in the formation of massive stars can be circumvented via an anisotropy of 
the thermal radiation field. Such an anisotropy naturally establishes with the forma- 
tion of a circumstellar disk. The required angular momentum transport within the 
disk can be provided by developing gravitational torques. Radiative Rayleigh-Taylor 
instabilities in the cavity regions - as previously suggested in the literature - are not 
required and are shown to be not occurring in the context of massive star formation. 



1 Introduction 

During their lifetime, massive stars exert a radiation pressure onto their surround- 
ings, which is higher than their gravitational attraction. How these massive stars can 
sustain accretion was an open question for decades. Semi-analytical (Kahn, 1974) 
and first radiation-hydrodynamical computations in spherical symmetry (Yorke and 
Kriigel, 1977) support the idea that the radiation pressure of a star retards the col- 
lapse of its proto-stellar core, leading to an upper stellar mass limit of < 40 M©. 

Based on the fact that in these ID simulations the collapse was reversed by the 
re-emitted radiation, Nakano (1989) inferred that the radiation pressure onto the 
accretion flow can be diminished by an anisotropic optical depth of the environment. 
Such an anisotropy can be produced by the formation of an accretion disk. 

The first numerical study aimed at verifying this idea was by Yorke and Sonnhal- 
ter (2002). In their simulations, the optical depth of the forming accretion disk 
turned out to be too low to provide a sufficient anisotropy of the radiation field. 
As a result, the accretion onto the proto-star stopped shortly after the disk forma- 
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tion, limiting the maximum final stellar mass to 42.9 M . In Krumholz et al. (2009), 
the authors claimed further feeding of the star and its circumstellar disk by a radia- 
tive Rayleigh-Taylor instability of the bipolar cavity regions. The most massive (and 
still accreting) proto-star formed in this simulation was 41 .5 M . 

All in all, no star much more massive than the ID radiation pressure barrier of 
M™ ax 40 M© was formed in numerical simulations up to this point. 



2 Methods 

In this section, we highlight two of our code specifics, namely the grid in spherical 
coordinates and the hybrid radiation transport method developed. For a complete 
description of the equations, methods, and numerical solvers please see Kuiper et al. 
(2010b) for the radiation transport, Kuiper et al. (2010a) for the hydrodynamics in 
axial symmetry, and Kuiper et al. (201 1) for the hydrodynamics in three dimensions. 

The grid of the computational domain is given in spherical coordinates with a 
logarithmically increasing resolution towards the center. The resolution of the inner 
region around the centrally forming star is 1 AU, decreasing towards the outer 
core regions proportional to the radius. This kind of grid covers the essential phe- 
nomena on the wide range of spatial scales from the outer core collapse (gravity 
dominated) over the torus formation (flattened rotating structure, still in-falling) in 
the transitional region down to the innermost disk formation (centrifugal balance). 

The hybrid radiation transport method is highly adapted to the problem under 
investigation: The stellar irradiation is computed via a very accurate frequency- 
dependent ray-tracing step, while the re-emission and thermal dust emission is com- 
puted within a fast gray Flux-Limited Diffusion (FLD) approach. 

Hence, the feedback onto the directly irradiated regions around the massive 
proto-star - the inner disk rim and the cavity walls - is computed with very high 
accuracy, including the wide range of different optical depths for the broad stellar 
spectrum. The ray-tracing of the stellar spectrum guarantees the resembling of the 
long range force of the irradiation, see Kuiper et al. (2012). 



3 Resolving the optical depth of the circumstellar accretion disk 

In our simulation series for varying sizes of the inner sink cell (Kuiper et al., 2010a), 
we find that it is essential to fully include and resolve the innermost part of the dusty 
accretion disk, which depicts the highest optical depth, to compute the anisotropy 
of the thermal radiation field correctly. 

Using the initial conditions of Yorke and Sonnhalter (2002), in which the size of 
the sink cell is several times larger than the dust sublimation radius, our simulations 
resemble their results of short disk accretion epochs. By restricting the sink cell to 
a size smaller than the dust sublimation radius, the thermal radiation field becomes 
strongly anisotropic and accretion sustains for several free-fall times. 
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4 Overcoming the radiation pressure barrier 

In simulation series for varying initial core masses M com in spherical as well as 
axial symmetry, we demonstrate, how the anisotropy of the thermal radiation field 
reduces the feedback onto the accretion flow, allowing the formation of the most 
massive stars known, see Kuiper et al. (2010a). In spherical symmetry, the accre- 
tion is stopped and reversed by radiative forces. Regardless of the initial core mass 
(A/core = 60 . . . 480 M Q ), the final mass of the forming star is limited to < 40 M . 
These results fully support the work by Kahn (1974) and Yorke and Kriigel (1977). 

Setting the initial pre-stellar core in slow solid body rotation yields the forma- 
tion of an accretion disk around the centrally forming massive star, which self- 
consistently leads to an anisotropy of the thermal radiation field. The optical depth 
of the disk remains at high values by further accretion from the large scale envelope. 
The bipolar cavity regions remain stable, and in turn lead to a mass loss of the proto- 
stellar core of roughly 50% of the initial core mass. In the shielded disk regions, the 
diminished radiative feedback onto the accretion flow enables the massive proto-star 
to increase its mass to M™ ax > 100 M . 

In fact, these are the first simulations, including the effect of radiation pressure 
feedback, which demonstrate a possibility to form stars up to the maximum value 
of the observed stellar mass spectrum. 

5 Angular momentum transport in massive accretion disks 

In the 2D simulations we have to rely on an a-viscosity model for the actual an- 
gular momentum transport. In Kuiper et al. (2011), we demonstrate in a three- 
dimensional simulation that the self-gravity of the forming massive accretion disk 
self-consistently leads to gravitational instabilities and the formation of spiral arms, 
which gravitational torques in turn yield an angular momentum transport. In con- 
trast to the smooth viscous disk accretion in axial symmetry, the 3D simulation 
shows episodic accretion. The accretion rate integrated over several episodic ac- 
cretion events is as high as (even slightly higher than) the accretion rates of the 
axially symmetric simulations in Kuiper et al. (2010a). Hence, we conclude that 
the accretion rate required to form massive stars is self-consistently arranged by the 
self-gravity of the forming massive accretion disk. 

6 On the radiative Rayleigh-Taylor instability 

In Krumholz et al. (2009) the authors claim that a 3D radiative Rayleigh-Taylor 
instability in the bipolar cavity shell is required to allow further feeding of the star- 
disk system beyond the radiation pressure barrier. Contrary, our simulations show 
the launching and expansion of stable outflow cavities, regardless of the dimension 
and for a large variety of initial conditions. Among others, both studies differ in the 
treatment of the direct stellar irradiation feedback. In Kuiper et al. (2012) it is shown 
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that if we cut down the radiation transport in our simulations to the gray FLD ap- 
proximation, the forming cavities undergo a radiative Rayleigh-Taylor instability as 
well. But if the stellar irradiation is computed with the much more sophisticated ray- 
tracing scheme, the cavities remain stable. This result is backed up in Kuiper et al. 
(2012) by analytical estimates of the radiative forces in the cavity shell, which are 
underestimated in the gray FLD approximation by up to two orders of magnitude. 

7 Summary 

In 45 published simulations regarding the radiation pressure feedback in the for- 
mation of massive stars, we have demonstrated a self-consistent and detailed de- 
scription of the formation of the most massive stars known. The radiation feedback 
in these simulations is computed via a sophisticated hybrid radiation transport ap- 
proach, highly adopted to the problem under investigation, see Kuiper et al. (2010b). 
On the one hand, the simulation series fully recover the previous ID radiation pres- 
sure barrier results of Kahn (1974) and Yorke and Kriigel (1977). On the other hand, 
the simulation series attest the numerical and physical improvements in our studies 
of the radiation pressure feedback in the multi-dimensional context with respect to 
previous numerical studies in the field, see Kuiper et al. (2010a, 2012). The radia- 
tion pressure feedback onto the disk accretion flow is diminished by an anisotropy 
of the radiation field due to the forming massive accretion disk, see Kuiper et al. 
(2010a). The required angular momentum transport will self-consistently be pro- 
vided by gravitational torques in the accretion disk, see Kuiper et al. (201 1). 

All in all, this demonstrates that the final masses of stars forming by accretion 
are not limited by the well-known radiation pressure barrier. 
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